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I.  Introduction 


So  far  the  existing  literature  on  the  multi  user  random  access  communication  systems  has  been 
dealing  with  a  homogeneous  population  of  users  [l]-[5].  There  are  many  practical  applications,  however, 
where  it  is  desired  that  some  packets  experience  shorter  delays  than  the  average  regular  packet  of  the 
system.  If  all  users  are  to  use  the  same  communication  system,  then  the  need  for  dividing  the  population 
of  users  into  two  classes  arises. 

There  are  cases  of  communication  systems  with  homogeneous  population  of  users  where,  at  specific 
known  time  periods,  the  input  traffic  to  the  channel  decreases  significantly  with  respect  to  the  nominal 
point  of  operation  of  the  system.  As  a  result,  the  average  packet  delay  decreases  but  the  utilization  of  the 
system  decreases  as  well.  Under  those  conditions,  we  can  improve  the  utilization  of  the  system  by  letting 
a  second  class  of  same  priority  users  have  access  to  the  system.  By  controlling  the  rate  of  the  input  traffic 
coming  from  the  second  class,  we  can  achieve  induced  average  packet  delays  for  both  classes  around  the 
nominal  point  of  the  original  class.  In  that  case,  the  same  algorithm  applies  to  both  classes  and  in  fact  we 
have  a  homogeneous  user  population.  A  second  option  is  to  adopt  an  algorithm  that  gives  priority  to  the 
packets  of  the  original  class.  In  that  case,  it  is  expected  that  if  the  induced  average  packet  delay  of  the 
original  (high  priority)  class  is  around  its  nominal  value,  then  the  induced  delays  of  the  second  class  will 
be  significantly  larger.  On  the  other  hand,  the  low  priority  packet  traffic,  that  induces  the  nominal 
average  packet  delay  for  the  high  priority  class,  is  expected  to  be  much  larger  than  in  the  previous  case  of 
the  equivalent  classes.  If  the  users  of  the  second  class  can  wait  for  the  occurence  of  the  low  traffic  time 
periods  of  the  original  system,  then  it  is  reasonable  to  assume  that  those  users  can  tolerate  an  additional 
delay  of  a  small  number  of  packet  lengths.  Thus,  by  using  a  system  with  users  with  different  priorities, 
we  can  greatly  increase  the  utilization  of  a  system  at  essentially  no  cost. 

In  a  mobile  user  environment  where  users  move  in  and  out  of  the  range  of  the  system,  or  move  from 
region  to  region,  fast  moving  users  may  need  to  experience  shorter  delays  than  the  regular  ones,  this  may 
be  necessary  to  make  packet  transmission  possible  while  the  user  is  still  inside  the  region.  Also,  users 
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that  arc  close  to  the  boundaries  of  a  region  and  are  going  to  move  outside  it,  should  experience  shorter 
delays. 

In  a  static  user  environment  there  are  also  cases  in  which  some  packets  have  high  priority  and 
should  reach  their  destination  faster  than  the  regular  ones.  High  priority  packets  can  be  those  which  are 
generated  by  high  priority  users  (e.g.  important  users,  or  users  that  can  pay  more  for  better  service),  or  can 
be  packets  that  are  generated  by  any  user  of  the  system  but  the  information  that  is  carried  is  characterized 
as  important  and  deserves  high  priority  in  its  transmission. 

An  important  measure  of  performance  of  a  communication  system  is  the  induced  average  packet 
delay.  In  some  environments,  there  may  exist  strict  constraints  on  the  delay  that  some  packets  can 
tolerate.  If  a  threshold  is  exceeded,  the  packet  is  considered  to  be  lost  and  the  average  number  of  those 
packets  can  be  a  measure  of  performance.  By  considering  that  those  special  packets  form  a  separate  class 
which  is  given  priority  by  the  system,  we  might  be  able  to  reduce  the  induced  delays  of  those  packets 
below  the  rejection  threshold  and  thus  greatly  improve  the  performance  of  the  system. 

In  the  next  two  sections  the  communication  system  and  the  suggested  algorithm  are  described.  In 
sections  IV  and  V  throughput  and  delay  analysis  arc  performed,  while  in  the  last  section  the  results  of  the 
analysis  are  shown  and  conclusions  are  drawn. 

II.  The  Communication  System 

We  consider  a  large  population  of  users  that  use  a  single  communication  channel.  We  assume  that 
users  which  for  some  reason  need  to  have  some  priority  over  the  rest  of  the  population,  form  the  high 
priority  class.  It  is  assumed  that  the  packet  traffic  generated  by  that  class  represents  only  a  small 
percentage  of  the  total  traffic  that  is  served  by  the  system.  In  other  words,  we  assume  that  the  packets  that 
need  special  service  arc  rare  and  this  is  a  realistic  assumption  at  least  for  the  environments  that  were 
described  above. 


The  input  traffic  to  the  channel  that  is  generated  by  each  class  of  users  is  assumed  to  be  Poisson 
distributed  with  intensities  Xj  and  X,  respectively;  the  Poisson  model  is  proved  to  be  an  appropriate  model 


for  the  cumulative  traffic  that  is  generated  by  a  large  population  of  bursty  users,  which  is  assumed  to  be 
the  case  in  the  system  under  consideration. 

Messages  are  assumed  to  be  packetized  and  of  fixed  length;  it  is  assumed  that  time  axis  is  slotted 
and  that  the  beginning  of  a  packet  transmission  coincides  with  the  beginning  of  a  slot. 

All  users  may  access  the  channel  as  long  as  they  have  a  packet  to  transmit;  the  first  transmission 
attempt  takes  place  at  the  beginning  of  the  first  time  slot  that  follows  the  packet  generation  instant. 
Because  of  the  freedom  that  the  users  enjoy  in  accessing  the  channel,  a  transmission  attempt  results  in 
either  a  successful  packet  transmission,  or  in  a  packet  collision  if  more  than  one  packet  transmissions 
were  attempted  in  the  same  time  slot.  Thus  it  becomes  obvious  that  an  algorithm  is  necessary  in  order  for 
the  conflicts  to  be  resolved  and  the  channel  to  remain  usable. 

It  is  assumed  that  all  users  that  have  a  packet  to  transmit  (and  only  these  users  need  to  do  that)  keep 
sensing  the  channel  and  are  capable  of  detecting  a  packet  collision;  that  is,  we  assume  that  a  binary 
feedback  information  is  available  to  all  active  users  before  the  end  of  the  current  slot,  revealing  whether 
the  slot  was  involved  in  a  packet  collision  (C)  or  not  (NC).  Channel  errors  are  not  taken  into 
consideration  and  packet  collision  is  the  only  event  that  results  in  unsuccessful  transmission. 

III.  Description  of  the  Algorithm 

The  first  time  transmission  policy  is  kept  the  same  for  both  classes  of  users;  it  is  simple  and  implies 
that  a  packet  is  transmitted  at  the  beginning  of  the  first  slot  following  the  packet  generation  instant.  It  is 
apparent  that  if  the  two  classes  are  to  experience  different  delays,  they  should  follow  different  steps  in  the 
collision  resolution  procedure. 

We  arc  going  to  use  a  simple  limited  sensing  collision  resolution  algorithm.  The  limited  sensing 
characteristic  is  apparently  important  for  a  mobile  user  environment  since  the  users  may  not  be  able  to 
know  the  history  of  the  channel  before  their  packet  generation  instant.  We  assume  that  the  state  of  a  user 
is  determined  by  the  content  of  a  counter  that  is  assigned  to  each  one  of  them;  this  counter  is  updated 
according  to  the  steps  of  the  algorithm  and  the  feedback  from  the  channel.  Users  whose  counter  content 
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at  the  beginning  of  a  time  slot  is  equal  to  one,  transmit  in  that  slot 

Let  Cjf(c*)  denote  the  counter  content  of  a  high  priority  (regular)  user,  at  the  beginning  of  the  i*  time 
slot.  Let  also  Fif  Fj£(C,NC),  denote  the  channel  feedback  information  just  before  the  end  of  the  i*  time 
slot.  The  steps  of  the  collision  resolution  algorithm  consist  of  the  following  counter  updating  procedures 
that  take  place  at  the  end  of  each  time  slot. 


(A)  IfFj  =  C  then 


(B)  If  F;  =  NC  then 


^  Cjf+1=l  with  probability  $ 
S*cLi~2  probability  l-$ 

^,c^,=2  with  probability  o 
with  probability  l-o 


c/  =  r  -»  ci,+1=rf2  ,  r£2  ,  je(s,f) 


c(  =  r  -»  c^,  =  r-1  ,  r£l ,  je(s,f) 

The  first  time  transmission  policy  can  also  be  described  by  using  the  concept  of  the  counter;  it 
simply  implies  that  a  new  user  sets  the  counter  equal  to  one  at  the  end  of  the  slot  in  which  its  packet 
arrival  took  place.  It  did  not  seem  to  us  reasonable  to  develop  different  first  time  transmission  policies  for 
the  two  classes  of  users.  It  would  probably  be  a  waste  of  the  channel  capacity  to  give  priority  to  rarely 
appearing  high  priority  packets,  before  it  becomes  known  that  a  collision  took  place.  If  a  conflict  occurs, 
then  the  collision  resolution  algorithm  offers  some  priority  to  the  high  priority  packets  that  were  involved 
in  the  conflict. 


From  the  description  of  the  algorithm  it  can  be  easily  obs  rved  that  the  system  is  of  continuous 
entry,  i.e.  new  users  enter  the  system  at  the  beginning  of  the  first  slot  that  follows  their  packet  arrival, 
unlike  what  happens  in  the  blocked  access  algorithms  |?|;  furthermore,  it  is  obvious  that  the  system  is  a 
last-come-first-served  one.  The  limited  sensing  characteristic  of  the  algorithm,  together  with  the  lack  of 
need  for  a  central  controller  to  coordinate  the  users,  increase  the  robustness  and  applicability  of  the 


system. 


IV.  Throughput  Analysis 


In  this  section  we  derive  bounds  on  the  stability  region  of  the  algorithm.  For  this  purpose,  we  use 
the  concept  of  the  session  and  develop  recursive  equations  to  describe  the  operation  of  the  system.  A 
session  is  defined  as  a  number  of  consecutive  slots,  as  it  is  explained  in  the  next  paragraph.  If  p  high 
priority  users  and  v  regular  ones  attempted  a  packet  transmission  in  the  first  slot  of  a  session,  then  the  pair 
(p,v)  determines  the  multiplicity  of  that  session. 

In  the  sequence,  we  define  a  session  of  multiplicity  (p,v),  p>0,  v>0  and  p+v>l,  by  using  the 
concept  of  a  virtual  stack  and  a  marker,  the  stack  is  assumed  to  have  infinite  number  of  cells.  We  assume 
that  the  system  starts  operating  at  time  t=0  and  that  the  marker  is  placed  in  cell  0.  The  first  slot  involved 
in  a  packet  collision  marks  the  beginning  of  a  session  of  multiplicity  (p,v),  if  p+v  packets  were  involved 
in  that  original  collision.  At  this  time  the  marker  is  placed  in  either  cell  3  or  cell  2,  depending  on  whether 
a  low  priority  user  was  involved  in  that  conflict,  or  not.  In  the  sequence,  the  market  moves  two  cells 
upwards  or  one  cell  downwards,  depending  on  whether  the  feedback  was  C  or  NC,  respectively.  The 
movement  of  the  marker  takes  place  at  the  end  of  a  slot.  The  slot  in  which  the  marker  moves  to  cell  0,  is 
the  last  slot  of  the  session.  The  first  slot  involved  in  a  collision  that  will  follow,  marks  the  beginning  of 
another  session  of  multiplicity  (p,v),  p>0,  v£0  and  p+v>l,  if  p  high  priority  and  v  low  priority  packets 
were  collided.  Sessions  of  multiplicity  (p,v),  p+v<l,  result  in  no  movement  of  the  marker  and  are 
defined  to  have  length  equal  to  one  time  slot.  It  should  be  noted  that  sessions  cannot  be  identified  by  the 
users  and  that  they  are  only  used  in  the  analysis  of  the  operation  of  the  system. 

From  the  previous  definitions  it  is  easy  to  conclude  that  the  multiplicities  of  the  sessions  arc 
independent  identically  distributed  random  variables  with  probability  density  function 
Pf,(p.v)  =  Pf(p)P,(v),  where  P,(F=p),  Ps(S=v)  arc  Poisson  density  functions  with  parameters  Xf  and  Xf 
respectively  ;  Xf  and  X,  arc  the  cumulative  input  traffic  rates  generated  by  users  with  or  without  priority 
respectively.  At  this  point  we  give  the  following  definition  for  the  stability  region  of  the  system. 

Definition: 
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If  for  an  input  traffic  pair  (Xf,  Xf),  the  expected  value  of  the  session  length  of  multiplicity  (|X,v)  is 
finite,  for  |X  and  v  finite,  then  we  say  that  the  operation  of  the  system  is  stable  and  the  pair  (Xf,  Xs)  belongs 
to  the  stability  region  of  the  system.  The  maximum  overall  sets  of  stable  points  (Xf,  Xs)  determines  the 
maximum  stable  throughput  region  and  is  denoted  by  S^. 

Let  t(p,v)  denote  the  length  of  a  session  of  multiplicity  (jx,v).  From  the  description  of  the  algorithm 
we  derive  the  following  recursive  equations. 

T0,0  =  X0,l  =  *1,0  =  1 

Vo=1+  X0,+Fi.s,  +  Vfc+Fa*  »  ( 1 ) 

X0,v  =  1  +  *F1(Sl  +  %0,-S,  +  *F„v^1+Sj.  V>2 

tH.v  =  1  +  *$,+F1(Si  +  V-$,+F„o,+S,  +  XF,,v-o,+S,'  1 .  1 

where  Fi(  S;  are  independent  Poisson  distributed  random  variables  with  parameters  Xf  and  Xs  respectively; 

<t)j,  Oj  are  independent  random  variables  that  follow  the  binomial  distribution  with  parameters  p,<!>  and  v,a 

respectively. 

Let  L^v  be  the  expected  value  of  the  length  of  a  session  of  multiplicity  (|i.,v).  By  considering  the 
expectations  of  both  sides  of  the  equations  in  (1),  we  obtain  the  infinite  dimensional  linear  system  of 
equations  with  respect  to  L^  v 

oo  oo 

L*v  =  Vv+IIa£V  Lkj  •  ^.v>0  (2) 

k=0j=0 

where  h^v=l  for  p£0,  v>0  and  a£’v  >  0,  for  all  nonnegative  p,v,k,j,  are  given  in  Appendix  A. 

Since  it  is  not  possible  for  the  above  system  to  be  solved,  we  will  try  to  find  bounds  on  L^  v.  The 
existence  of  an  upper  bound  on  L^v,  that  is  finite  for  p  and  v  finite,  will  provide  a  lower  bound  on  the 
maximum  stable  throughput,  according  to  the  previous  definition.  We  found  that  the  quantity 

i  u  -i  u  =LU  =1 

MJ.O  ~  Lo,i  Lt,o  1 

L“v  =  ap  +  pv  +  y,  l<p+v<°°  (3) 


satisfies  the  inequality 


(4) 


Vv+IZa^VL;iV<L;v  l<n+v<=o 
k=0j=0 

for  a  =  a(Xf,Xs,<p,o),  p  =  (3(Xf,Xs,<j>,c)  and  y=y(Xf,Xs,(j),a)  properly  chosen.  To  find  proper  values  for 
a,  P  and  y,  we  substituted  (3)  into  (4)  and  three  categories  of  infinite  number  of  inequalities  were 
obtained.  By  properly  choosing  yand  p  to  satisfy  some  initial  inequalities  and  by  using  numerical  search 
techniques  and/or  the  limiting  inequalities  as  p— and/or  we  found  a  proper  value  for  a,  so  that 
all  three  categories  of  inequalities  be  satisfied. 

It  is  trivial  to  enumerate  the  unknown  quantities  L^v  that  appear  in  (2),  by  finding  a  mapping  rule 
from  the  set  Q  =  ((p,v),  peZ^,  veZj)  onto  the  set  of  nonnegativc  integers  Z^.  Then,  the  system  in  (2)  can 
be  written  as 

oo 

Li  =  hi  +  Z,  ajLj  (5) 

j=o 

where  ieZj  corresponds  to  a  specific  pair  (p,v).  The  references  that  are  given  throughout  the  analysis  are 
related  to  the  solutions  of  systems  of  the  form  (5).  Since,  as  we  have  just  noticed,  systems  (5)  and  (2)  arc 
equivalent,  the  results  that  appear  in  the  references  concerning  the  solutions  of  (5)  extend  naturally  to  the 
solutions  of  (2). 

From  the  existence  of  the  L“v  and  since  h^v  >  0  and  a^v  >  0.  for  all  nonnegativc  p,v,k,j,  it  is 
implied  ([7],(  1  ])  that  the  infinite  dimensional  linear  system  of  equations  (2)  has  a  unique  nonnegativc 
solution  that  satisfies 


0  <  L^v  <  L“v  =  ap  +  Pv  +  y. 

In  Fig.  1 ,  a  lower  bound,  S^ax,  on  SmiX  is  plotted;  the  stable  region  includes  all  points  (\,\)  for  which  a 
bound,  as  in  (3)  that  satisfies  (4),  was  possible  to  be  obtained. 


If  only  the  high  priority  class  is  using  the  communication  channel  then  the  system  will  be  able  to 
serve  larger  traffic  (Xf  <  0.357)  than  in  the  case  in  which  only  low  priority  users  were  served  (Xs  <  0.32). 
This  is  not  surprising  since  the  procedures  of  the  algorithm  arc  the  same  for  both  classes,  with  the  only 
significant  difference  that  there  is  always  a  waste  of  the  first  slot  after  a  collision  among  low  priority 


users.  As  a  coasequence,  there  is  a  reduction  in  the  maximum  stable  throughput  when  only  the  class  of 
the  low  priority  users  is  served. 

An  upper  bound,  on  the  maximum  stable  throughput  can  also  be  obtained.  This  bound  is 
given  by  the  set  of  all  pairs  (Xf,  A.s)  that  provide  a  nonnegative  solution  to  a  truncated  version  of  the 
system  in  (2)  [8],  [4],  By  using  a  large  number  of  equations  we  generally  obtain  a  tight  upper  bound. 
The  upper  bound  that  was  obtained  for  p<5,  v<10  is  plotted  in  Fig.  1. 


V.  Delay  Analysis 


In  this  section,  upper  and  lower  bounds  on  the  average  packet  delay  are  derived.  It  turns  out  that,  on 
the  average,  high  priority  packets  experience  shorter  delays  the  regular  ones. 


The  existence  of  renewal  slots,  under  stable  operation  of  the  system,  that  marie  the  beginning  of 
statistically  identical  sessions,  implies  that  the  operation  of  the  system  can  be  described  by  a  regneraiive 
process.  Under  these  conditions,  we  can  draw  conclusions  about  the  limiting  behavior  of  the  system  by 
manipulating  quantities  that  are  defined  on  a  session  [9],  [10],  The  application  of  regenerative  theory 
procedures  to  the  delay  analysis  of  random  access  algorithms,  appears  in  [12],  [4],  The  same  results  can 
be  obtained  by  using  directly  the  strong  law  of  large  numbers  [11],  [7], 


Let  Wf  and  Ws  be  the  mean  cumulative  delay  of  all  high  priority  or  regular  packets  respectively, 
that  arrive  in  a  single  session;  the  time  interval  between  a  packet  arrival  instant  and  the  beginning  of  the 
time  slot  that  follows  the  packet  arrival,  is  not  included  in  Wf  or  Ws.  If  Df  and  Ds  denote  the  average 
delay  of  a  high  priority  or  a  regular  packet  respectively,  then  from  the  discussion  in  the  previous 
paragraph,  we  have 


Wf'  W,u 

0.5  + - <  Df  <  0.5  + - -  (6a) 

XfC  \Ll 


w '  w“ 

0.5  + - <  Ds  <  0.5  + - -  (6b) 

XSLU  XSL 

where  Wf,  w],  L1  and  Wfu,  W“,  Lu  denote  lower  and  upper  bounds  on  the  corresponding  quantities;  L  is 


wwswwnii  wwwwwwwfwwwnww  im  rmmn 


the  average  sessions  length  and  .5  is  the  mean  packet  delay  until  the  beginning  of  the  first  time  slot  that 
follows  the  packet  arrival. 


A  lower  bound  on  L  can  be  obtained  by  solving  a  truncated  version  of  the  system  in  (2).  Let 
be  the  solution  of  the  finite  system  in  which  only  L^v  for  p<M,  v^N  are  considered.  For  (>^,  it 

is  guaranteed  that  <  L^v  ([8],  [7])  and  thus,  is  a  lower  bound  on  L^v.  By  assuming 

L^vN  =  0,  |i>Morv>N,  and  considering  the  expected  value  of  L^1/4  with  respect  to  (p,v)  we  obtain  a 
lower  bound,  L1,  on  the  mean  session  length. 


The  quantity  that  appears  in  (3)  can  serve  as  an  upper  bound  on  L^v.  This  bound  is  arbitrary  and  it 
r  m.n  . 


is  generally  loose.  If  L^v  is  the  solution  of  the  finite  system  of  equations 


where 


f  M  N  _  r  M.N  _  r  mjm  _ 

L0.0  “  L0,l  ~  LI.O  ~  1 

CN  =  +  X Z  akjV]-kj’N  0<p<M,  0<v<N,  1  <p+v 

k=OH> 


*CN  =  V  +  I  I  «£vLi 


(7a  i 

(7h) 


k=M+1  j=N+1 
i 

kj’  ~  &•  \~,tt - “|i,v  - -rr-  — - - 

L^y  14],  (8).  An  upper  bound  on  L  is  thus  obtained  by  considering  the  expected  value  of  L^v  with  respect 
to  (p.,v);  where 


and  L^:,  k>M  or  j>N  are  given  by  (3),  then  L^f*  is  an  upper  bound  on  L  that  is  generally  tighter  than 


L“v  =  C"’  °-^M  -  °-v-N  (8a) 

L^v  =  L£v>  otherwise.  (8b) 

In  the  sequence,  we  derive  bounds  on  Wf  and  V.'s.  Let  co^v  and  0)^v  be  the  cumulative  delays  of  the 
high  priority  and  the  regular  packets,  respectively,  that  arrive  during  a  session  of  multiplicity  (p,v).  It  is 
easy  to  observe  from  the  description  of  the  algorithm  that  0)^v  and  G)*v  satisfy  the  following  recursive 
equations. 

^0,0  =  ®0,1  =  ®  ’  0}|,0=^  (9a) 
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vrar-l 

yVvi 

*  »  *  »  » 


,\.v 
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“£0  =  ^  +  <4^  +  (M--4>i>tr  ^,.s,  +  “U+F,s,  •  ^2  (9b) 

“av=“F1.SI  +  «4J.ol+Sa  +  4J,v-a1+S1  •  VS2  (9c) 

<v  =  H  +  <+fiA  +  (HH>,^i+FiiSi  +  co^i+Fj0i+S5  +  ^  .  ^l.v>l  (9d) 

and 

®o,o  =  ®1,0  =  0  >  %1 =  1  (9e) 

<0  =  <+F,.S,  +  “^1+F2^  .  (9f) 

®0,v  =  V  +  <1(Si  +  VTFl>Sl  +  <,o1+Sj  +  +  <>v^1+s,  .  v>2  (9g) 


<V  =  V  +  <+FliSl  +  vt0i+fiiSi  +  (V-O^T^f^s,  +  “(U.^o.+S,  +  <.v^1+S,  .  ^1,V>1  (9h) 

where  all  variables  are  as  defined  in  (1).  By  considering  the  expectations  of  the  above  equations  we 

obtain  the  following  infinite  dimensional  linear  systems  of  equations 


A 

o 

II 

gl 

II 

o 

S: 

r- 

o 

II 

»— * 

(10a) 

oo  oo 

=  •  I«*>1 

(10b) 

k=OpO 

W‘o  =  W‘o  =  0,  wj,-i 

(10c) 

oo  oo 

K.-si.+  EZaffXj  •  ‘‘+v>1 

(lOd) 

k=0j=0 

where  g£v,  g^v  are  given  in  Appendix  B  and  a£]v,  for  p.,v,kj  nonnegative,  are  the 

same  as  in  (2)  and 

appear  in  Appendix  A. 

By  following  procedures  similar  to  those  that  were  used  to  obtain  bounds  on  L,  we  derive  upper  and 
lower  bounds  of  Wf  and  Wr  It  turned  out  that  a  universal  initial  upper  bound  on  W^v  and  W^v  was  hard 
to  obtain.  Thus,  we  had  to  divide  into  several  regions  and  derive  bounds  valid  for  (Xf,Xs)  in  a 


specific  region.  At  that  point,  we  made  the  assumption  that  the  high  priority  traffic  is  much  lower  than 
the  regular  one  and  assumed  X,-  to  correspond  to  less  than  20%  of  the  total  traffic.  This  assumption 
seemed  to  us  to  be  realistic  and  bounds  were  obtained  for  the  operation  region  of  the  system 


Sop  H  ( W  :  0^^0.065, 


Vnax'* 


where  X.  m„(Xf)  can  be  obtained  from  Fig.  1.  The  following  upper  bounds  on  W^v  and  W^v  satisfying 
inequalities  similar  to  those  in  (4),  were  obtained. 

w£v  =  S(J-2  +  ev2  +  0|X  +  L  “ v  ,  O^SO.065 

w£“  =  8p2  +  ev  +  L  “>v  ,  0<k<k*(kf) ,  Oi^SO.065 
wi“  =  5p2  +  ev2  +  0v  +  L“v ,  X**(^X^\mJXf) ,  0<Xf<0.065 
W  =  5p2  +  ev2  +  L “ v  ,  0<Xs<Xs**(Xf) ,  0<Xf<0.065 


dla) 

(Hb) 

(He) 

(Hd) 


where  Xi*(Xf)  and  Xg**(Xf)  are  some  values  of  X,  less  than  Xi  rrax(Xf)  and  L  “ v  is  given  by  (3). 

The  existence  of  the  above  upper  bounds  guarantees  that  the  finite  dimensional  linear  systems  of 
equations 


M  N 


W|^'N  =  g)I,v+ZZafiVwffiN  •  0£^M,  (KvSN 


k-Oj-0 
M  N 


(12a) 

(12b) 


W1I^  =  g;,v+ZZagvW 0^v5N 

k-0j-0 

have  a  unique  nonnegative  solution  that  is  a  lower  bound  on  W^v  and  W*v  respectively,  for 
0<p£M,  0£v<N.  By  using  zero  as  a  lower  bound  on  v  and  W^v  for  p  >  M  or  v  >  N,  and  considering 
the  expectations  with  respect  to  (p,v) ,  we  obtain  a  lower  bound  on  Wf  and  Wg  respectively. 


A  tighter  upper  bound  on  Wf  and  Wg  can  be  obtained,  as  explained  earlier  in  this  section,  by  solving 
the  finite  dimensional  linear  systems  of  equations 


-  f,M,N  _  -  f,M,N  _  _  ii/  _  • 

w0,0  —  w0,l  -  u  •  w|,0  -  * 


M  N 


=  Giv  +  II  aftv  W™  .  0<p<M,  0<v<N,  p+v>l 
k-Oj-O 


(13a) 

(13b) 


and 


W, 


J.M.N 


0.0 


=  w 


0  -o  ,  W01  -  I 


(13c) 
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(13d) 


M  N 

Wj**  =  G ‘  v  +  £ X  a£jv  wjf",  0<^M,  0<v£N,  p+v>l 

k-Oj-O 

where  G^v  and  G^v  are  given  in  Appendix  C.  An  upper  bound  on  Wf  and  W,  is  then  obtained  by 
considering  the  expected  value  of  w£“  and  w£“  with  respect  to  (p,v),  where 

W^  =  W^.  0^v<N 

w£“  =  w£“  ,  otherwise 
and 

0<m£M,  0Sv<N 
Wjj-Wj;  ,  otherwise 

The  bounds  on  L,  Wf,  Wf  that  were  obtained  in  this  section  for  some  values  of  (Xf,Xt)eSop  are  shown  in 
table  1.  By  substituting  those  bounds  in  (6),  we  calculated  tight  upper  and  lower  bounds  on  Dr  and  Ds; 
these  bounds  appear  also  in  table  1  for  some  values  of  (Xf,  X,)eSop. 

VI.  Results  and  Conclusions 

The  algorithm  that  we  developed  and  analyzed  is  supposed  to  operate  in  an  environment  where  two 
classes  of  users  with  different  priorities  are  accommodated.  An  algorithm  for  a  homogeneous  user 
population  that  would  work  in  a  similar  way  and  use  binary  feedback  information  and  simple  splitting 
after  a  collision,  has  been  found  to  achieve  a  maximum  stable  throughput  of  -  .36  [13].  The  algorithm 
that  we  suggest  for  a  non-homogeneous  population  achieves  total  throughput,  at  least,  between  .320  - 
.357  depending  on  the  contribution  of  the  two  classes  to  the  total  input  traffic. 

In  Fig.  2,  Fig.  3  and  Fig.  4,  plots  of  the  bounds  on  Df  and  Ds  versus  X,,  for  XpO.Ol,  X,=0.03  and 
Xf=0.065  respectively,  are  shown.  These  values  of  Xf  correspond  to  an  input  traffic  coming  from  the  high 
priority  class  equal  to  -  3%,  -  10%  and  ~  20%  of  the  total  traffic  that  can  be  served  by  the  system.  From 
the  plots  it  can  be  observed  that  the  high  priority  packets  experience  shorter  delays  than  the  packets  of  the 
other  class;  the  difference  is  essential  for  X,>.5X,  If  the  nominal  point  of  operation  of  the  system  is 
set  around  X,=.9Xi  nMX,  then  the  average  high  priority  packet  delay  is  less  than  half  the  one  of  the  other 


In  table  1,  the  delay  results  of  the  suggested  algorithm  are  compared  with  the  delay,  D*.  that  the 
homogeneous  class  equivalent  algorithm  (as  described  above),  induces  [13].  Again  we  can  observe  that 
always  Dj<D*  and  particularly  Df<.5D*  around  the  nominal  point,  the  latter  being  defined  as  before. 

Since  privileged  service  is  offered  to  some  users,  there  has  to  be  a  price  that  the  rest  of  the 
population  must  pay.  The  first  consequence  is  the  small  reduction  in  the  total  throughput,  as  mentioned 
before.  The  other  penalty  is  the  increased  average  low  priority  packet  delay  compared  with  the  one  that 
the  homogeneous  population  equivalent  algorithm  induces.  From  table  1  we  can  see  that,  indeed,  DS>D*. 
as  it  was  expected.  The  increase  in  Ds  is  far  from  catastrophic  and  it  is  realistic  to  consider  that  it  is 
possible  for  a  system  to  tolerate  these  delay  increases  for  the  low  priority  class,  especially  if  strict 
limitations  exist  for  the  high  priority  users. 

As  an  example,  consider  the  communication  system  described  in  the  second  paragraph  of  the 
Introduction.  Assume  that  the  input  traffic  of  the  original  class  at  the  nominal  operating  point  is  .25 
packets/packet  length  and  thus  the  (desired)  induced  average  packet  delay  is  5.5  -  6.0  packet  lengths  (last 
column  of  table  1).  Assume  that  at  night,  the  input  rate  falls  to  0.065  packets/packet  length.  At  that  time, 
a  second  class  of  users  is  given  permission  to  use  the  channel.  If  the  induced  average  packet  delay  of  the 
original  class  has  to  be  at  most  =  6.00  packet  lengths,  then  depending  on  the  case  we  observe  the 
following:  (a)  If  the  second  class  has  the  same  priority  as  the  original,  then  the  additional  input  traffic 
rate  that  can  be  accommodated  by  the  system  is  0.185  packets/packet  length,  (b)  If  the  second  class  has 
low  priority,  then  the  additional  input  traffic  rate  becomes  .25  packcts/packct  lengths  (table  1).  Thus, 
there  is  an  increase  by  =  35%  of  the  additional  traffic  that  can  be  accommodated,  if  the  population  of 
users  is  divided  into  two  classecs  with  different  priorities.  The  increase  in  the  average  packet  delay  of  the 
low  priority  users  is  rather  negligible  compared  to  a  realistic  waiting  time  until  these  users  are  given 
permission  to  access  the  channel. 


.1 
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Figure  1.  Upper,  SU  ,  and  lower, ,  bounds 

max  ’  max 

maximum  stable  throughput;  \  and  \ 
in  packets/packet  length. 
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Figure  3.  Average  packet  delay  of  the  high,  D^,  and  the 
low,  Dg,  priority  classes  (in  packet  lengths), 
versus  the  total  input  traffic  rate,  (in 

packets/packet  length).  For  *  .03. 
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Appendix  A 


The  coefficients  of  the  linear  system  of  equations  that  appear  in 
(2)  are  given  by  the  following  expressions. 


(a)  For  y=0,  v=0  or  u=0,  v=l  or  p=l,  v=0, 


“kj  =  ° 


for  all  k  _>  0,  j  >_  0. 


(6)  For  y=0,  v>2. 


(i)  for  0  _<  k  <  °°,  0  £  j  _<  v 


Pf<k>Ps«>  +  Pf(k)  J0  »J1)ps<3-1>  +  Pf<»  J0  bv(v-l)Ps(j-l) 


(ii)  for  0  £  k  <  <»,  v_<j  < 


V  V 

Pf(k)p  (j)  +  P  (k)  l  b  (v-l)P  (j-v+1)  +  pf(k)  l  b  (1)P  (j-v+1) 
J  S  r  1=0  V  1=0  v  S 


(y)  For  p  >2,  v  =0, 


(i)  for  O^k^p,  0j<j  <  °° 


^-VJ)  l  b  (k-i)P  (i)  +  P  (j)  l  b  (u-k+i)Pf(i) 
i=0  U  r  s  i=0  U  f 


(ii)  for  p  k  <  0  <  i  «■ 


u  p 

-£•{«  p8(j)  l  b  (i-i)Pf (k-p+i)  +  P  (j)  y  b  (i)P  (k-p+i) 
Kd  i=0  U  *  S  i=0  -  f 


(6)  For  y  >  1,  v  >  1, 


(i)  for  0_^k^u,  0  _<  j  _<  v 


k  k  i 

ak’i  =  ps(j)  l  b  (k-i)P  (i)  +  l  b  (u-k+i)P  (i)  J  b  (1)P  (j-1)  + 
k*J  s  1.-0  w  t  i=0  U  1  1=0  V  S 

j 

+  P  (k)  l  b  (u-l)P  (j-1) 

1  i  V  S 


(ii)  for  u^k<°°,  0  j  v 


u  u  j 

\’V.  =  pc(j)  J  b  (u-i ) Pf (k-u+i )  +  7  b  (i)P  (k-u+i)  >  b  (1)P  (j-1)  + 

3  i=0  u  1  i=0  1  1=0  v  S 


i 

+  P  (k)  l  b  (v-l)P  (j-1) 
r  1=0  u 


(Hi)  for  0  k  u,  v  <_  j  <  » 


k  k  v 

au’-  *  pc(j)  [  b  (k-i)P.(i)  +  l  b  (u-k+i)P,  (i)  7  b  (v-l)P  (j-v+1)  + 

k’J  i=0  11  f  i-0  U  1=0  '  s 

V 

+  pf(k)  l  b  (1)P  (j-v+1) 

1=0  v  s 


(iv)  for  u_^k<®,  v  _<  j 


U  U  V 

=  P  (j)  I  b  (y-i)P  (k-u+i)  +  7  b  (i)P  (k-u+i)  y  b  (v-l)P  (j-v+1)  + 

k*j  s  i=0  U  1  i  =  0  U  1=0  v'  S 


+  P  (k)  V  b  (1)P  (j-v+1 ) 
i=0 
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b  (1) 
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l!(v-l)! 


U,v 


The  following  summations  involving  the  coefficients  a  ’ .  were  of 

K  i  3 


wide  use  in  the  analysis  of  the  algorithm. 


oo  oo 


TO  -X) 


l  l  -“’v  -  ’  .  •  r 


=  3  (y>l,  v>l),  )  >  a,  \  =  2  (y>2) 

k=0  j=0  K,J  k=0  1=0  ,J 


OO  OO 


»JO  JO 


l  l  =  3  ('--2)  ’  i  l  ak’>  =  3vf 

k=0  j=0  k=0  j  =0 


+  u  (-..  I,  >n 


CO  oo 


OO  JC 


l  l  =  3Xc  +  V  (u-1,  V>1)  ,  l  7  aj^k  =  2-  +  y  (y>2) 


-.0, 


k=0  j=0 


k,jJ  s 


k=0  j  =0 


k,j  f 


oo  oo 
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=0  j=0  s  k=0  j=0  ,J 


jo  --o 
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OO  oo 


l  l  =  3Xs 

k=0  j=0 


'"40  'JO 


+  V  (V22)  ,  l  >  a£’V  =  2 ( '  + 

k=0  j=0 


>;>  (ve2) 


00  oo 


l  l  \’°.k2  =  2  ( >  +  X 1)  +  U2(l  -  23(1-3))  +  ii  ( 2  X  f  +  23(1-3))  (;^2) 

k=0  j=0  K,J  r  1 


2,  .  .  2, 


OO  oo 


l  l  -°’v'-2 


3,  7k  =  3  ( •  ,  +  K)  (v->2) 


k=0  j=0 
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For  the  derivation  of  the  upper  bounds  on  W  and  W  that  are  given 

u,v  u,v 


by  (11),  we  calculated  the  quantities  ^  and  g  ^  by  using  the  upper  bound 


(3)  in  the  place  of  ^  uPPer  bounds  (11)  are  also  upper 


bounds  on  the  solutions  of  the  system  in  (10),  where  the  exact  values  of 


f  s 

L  are  used  in  the  derivation  of  g  and  g  [7],  [8],  [14].  By 
p,v  °y,v  u,v 


using  the  upper  bound  in  (3)  we  found  the  following  expressions. 


(B  )=-->gf  «  a<Kl-4)p2  +  U  +  «  +  (1- P)  (otX  +  sx  +  y)  -  n>  - 

l  U ,  v  is 


-  a(b  (l)Pf (0)P  (0)  +  b  (0)Pf(l)P  (0))  -  £b  (0)Pf (0)Po(l) 
Uts  uts  ufs 


ot<f» (1— d>)  ]b  -  (a+B+l)b  (1)P,(0)P  (0) 

ufs 


where  fi  =  b  (0)Pf(0)P  (0)  +  b  (l)Pf(0)P  (0)  +  b  (0)P,(1)P  (0)  +  b  (0)P,(0)P  (0) 
Uts  Ufs  uts  uts 


(B2)  g^  ^  =  8a(l-a)u2  +  {(2-o)(aX^  +3Xg  +  y)  -  80(1-?)  + 


+  (l-y)(bj(0)Pf(0)Ps(0)  +  Pf (0)Pg (0) )  +  (1-S-Y)(b  (l)Pf(0)Ps(0)  + 


+  bv(0)Pf(0)Ps(l)  +  Pf  (0)PS(1))  +  (l-'i-Y)  (bv(0)Pf  (l)Pg(0)  + 


+  Pf(l)Pg(0))  +  l}v  +  (6-Y-l)bv(l)Pf(0)Ps(0) 


(B,)  _>  gf  *  B *o (l-o) v  +■  ai($+(l-:)  (1-  j  +  l(2-’)(  «‘  +  V  +  ,)  - 

j  p  $  v  r  s 


-  Bo(l-o)  +  ( 1— y  )  (b  (0)P  (0)b  ((.)P,(0)  +  b  (0)P.(0)P  (0))  + 

v  s  u  t  ufs 


+  (1-oc-y)  (b  (0)P  (0)  [P  (l)b  (u)  +  P,(0)b  (u-l)  ]  +  b  (0)P,(1)P  (0)  + 
v  s  i  u  r  u  u  i  s 


_  1*,"  f 


»vr 

»V,  i*’ 


I* 


& 
kV.V 


m 


p 

»)v' 


+  b(l)Pf(0)P  (0))  +  (1-B-Y)(b  (1)P  (0)b  (u)P, (0)  + 

Mis  V  S  M  I 


+  b  (0)Pf(0)P  (1)  +  b  (0)P  (l)b  (u)P,(0))  +  l}v  + 

P  1  S  V  5  y  X 


+  bv(l)Ps(0)by(uOPf(0)(i3+Y-l) 


f  g 

For  the  derivation  of  the  lower  bound  on  W  and  W  that  is  given  by  the 

u  ,v  p ,  v 


solution  of  (12)  we  simply  substituted  in  (B^) ,  (B2)  and  (B^) ,  the  lower 


bounds  on  L 

U ,  V 


Appendix  C 
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The  constants  ^  and  ^  of  the  systems  in  (13)  are  given  by 


Gl,0  "  1 


GJj  v  3  I  l  »  WO.  v>0. 

k=M+l  j=N+l  K,J 


u,v  k 


l  l  ak’i  wk’j  +  w  +  e^(m-*i)ta  +f  s 
=m  j-N+i  k,j  1  u  V  1 


}  ,  U>1  ,  v>0 


l  l  W?’“  .  U>0  ,  v=0 

k=M+l  j-N+1  K,J 


l  l  ak’X’U+  v+E(vtf  s  +  ( 
■M+i  i=n+i  r  1 


,  4.0  K  ^=0  .  V>1 

2’°1  S2 


Gs  -  V  T  a,U,v  Ws,u  +  V  4-  E  {  vt  _  +  (v-1.)  x  } 

U,V  k-W-1  j-N+1  k’^  k’^  VF1,S1  P-«1+F2,o1+S2  , 
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kri  ■ 
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m 
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in  the  place 


By  substituting  the  values  for  V/f ’ y  and  wf’f  from  (11)  and  L 

k,j  k,j  u,v 

from  (8),  we  finally  obtain 


Gi,o  "  1 


l  l  ?  !  ak*j  Wk’j 

k=0  j=0  k’J  k’j  k=0  j-0  k,J  k,j 


Gu  v  *  I  I  ak’j  wk’j  -II  aJ*’j  +  v  +  vi (1— <(>)  [o<(>(y-l)  + 
y’  k=0  j-0  ,;I  k=0  j-0  k,J  k,j 

+  ctA  +  BA  +  y]  -  £  l  P  (F  )P  (S  )b  ((j>  )  ta(<{,  +F  )+6S  +v-Ly  „  ]  (u-4>, 

VFi-M  1  VF1’S1  1 


,  u=0  ,  v>0 
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$ 

vv' 


!«s* 

■ft' 

>’;V 


I 


p>0  ,  v>0 


Go,i  -  1 


00  00  M  N 

GS  =  I  I  ay’y  wf»“  -  y  y  ay’y  Ws’u 

M,V  k=0  j-0  k’J  k,j  k=0  j-0  k,j  k,i 


U>0  ,  v=0 


M  N 


:',V  ”  l  l  l  l  wj’“  +  *  ♦  vl«xt+6x  *,] 

k-0  j-0  ,J  K,;l  k-0  j-0  k'J  1  8 


+  v(l-o) [aA^+8o(v-l) 


+  e»8  +  vl  +  I  J  Pf(F1)Ps(S1ML" 

r .  <M  11 


“Fr  Pf<F2)P,<S2)*-v«’1)(''-»l)lkF2,l-1+S2 


Oj+S^N 


aF2  -  B(o1+S2)~y]  ,  u=0  ,  v>l 


m 

y 


■S' 

»*»  -N 
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oo 


l  l 

k-0  j=0 
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y,v 

k,j 


w; 


,s  ,u 
k,j 


M  N 

I  I  Wf  +  v  +  v[ay<HaX  +BX  +y] 

k=0  j=0  f  s 


v(l-a)  [ay-ay<)>+aX  +ga(v-l)+BX  +y]  +  £  £  P  (F  )P  (Sn)b  (<J>  ) 


4»1+F1_<M 

S^N 


ill  p-(j!1+F„<M 


12- 
*!+S2lN 


[LU 


U~4'l+F2’al+S2  « Cp-<J>1+F2)  -  6(o1+S2)-y]  ,  u>0  ,  v>0 


The  calculation  of  the  infinite  summations  is  carried  out  by  using  the 
infinite  summations  that  appear  in  Appendix  A. 
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